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ABSTRACT 

A code, Epic5, has been developed which computes, in the two-dimensional case, the 
initially circular orbits of guiding centra in an arbitrary axisymmetric potential with an arbi- 
trary, weak perturbing potential in solid body rotation. This perturbing potential is given by 
its Fourier expansion. The analytic solution solves the linear epicyclic approximation of the 
equations of motion. To simulate the motion of interstellar matter and to damp the Lindblad 
resonances, we have in these equations introduced a friction which is proportional to the devi- 
ation from circular velocity. The corotation resonance is also damped by a special parameter. 
The program produces, in just a few seconds, orbital and density maps, as well as line of sight 
velocity maps for a chosen orientation of the galaxy. 

We test Epic5 by comparing its results with previous simulations and observations from 
the literature, which gives satisfactory agreement. The aim is that this program should be 
a useful complement to elaborate numerical simulations. Particularly so are its abilities to 
quickly explore the parameter space, to construct artificial galaxies, and to isolate various 
single agents important for developing structure of interstellar matter in disc galaxies. 
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1 INTRODUCTION 



Since the emergence of high spatial resolution two-dimensional 
velocity field for galaxies, the interpretation of observed velocity 
fields has become one of the key concerns of the extragalactic sci- 



entific community (e.g., Boulesteix et al.|1987| Adam et al. 


1989 


|Bacon et al.|1992||Bacon et al.|2001| Allington-Smith et al. 


1997 


Hernandez et al. 2003 1. Of particular importance is the class of 



galaxies which host bars and/or spiral arms, as they are estimated to 



be the most abundant types of galaxies in the local Universe jLove-| 
|day|1996) . 

The theoretical foundation for understanding the kinematics 
and dynamics of spiral galaxies has been worked out over the last 
century. Bertil Lindblad (1963 1964) conceived a picture of cir- 
culation of stars between the spiral arms of a quasi-steady rotat- 
ing spiral potential. Simultaneously, [Lin & Shu] < |1964| > presented 
their density wave theory approaching the problem from a differ- 
ent theoretical point of view. Subsequently, |Shu et ah] ( |1973| > and 
|Roberts et al.| ( |1979| > derived solutions for the circulation of the 
interstellar medium through such a spiral potential and predicted 
large scale galactic shocks along the spiral arms. These processes 
make bars and/or spiral arms potential actors to redistribute angu- 
lar momentum which will lead to dramatic effects such as intense 
star formation, build up of bulges, or onset of nuclear activity. In 
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turn, the host galaxies may undergo dramatic secular evolution on 
time-scales compared to the galaxy dynamical time (see Kormendy 
|& Ke nnicutt 2004]for a comprehensive discussion). 

Beginning with Holmberg ( 1941 1 in the pre-computer era and 
in the 1950s with electronic computers (e.g. P.O. Lindblad 19601, a 
large number of simulations of development and evolution of struc- 
ture in galaxies have been made with a variety of computer codes. 
These studies have often aimed at establishing the details of the 
flow of gas in and around bars (eg. Athanassoula 1992}, as well 
as the mechanisms which trigger starbursts and nuclear activity in 
the centers of galaxies (e.g., Shlosman et al. 1989). However, they 
often involve extensive computational power leading to relatively 
long computing times. 

Several successful attempts have been carried out in order to 
deliver analytic solutions for observable parameters for some of the 
effects involved in the above questions (Sakhibov & Smirnov 1989; 
Canzian|[T993l |Wa&jT994l |Lindblad & Lindblad|[T994"l |Schoen-| 
makers et al.|1997||Wong et al.|2004||Fathi et al.|2005]|Byrd et"aI1 
2006fl van de Ven & Fathi|2010} , though they have all exclusively 
treated one specific morphological feature in each analytic model. 
While a linear analysis can be a guide for a physical understanding 
of complicated numerical results from simulations, in these analytic 
studies, bars and spiral arms have been accounted for separately, 
and interpretations have been based on marginalizing the effect of 
the other features. A notable difference between these results and 
numerical simulations is that the stepwise integration in a detailed 
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simulation requires appreciable computing times, while for analyt- 
ical calculations the computing time is negligible. 

In this paper, we present an analytic solution within the 
epicyclic approximation, in which we introduce an arbitrary grav- 
itational potential and a damping coefficient for adequate appear- 
ance of the corotation resonance. Our solution has been coded as 
a development of the program Epic (Lindblad & Lindblad jl994) 
which was based on the first order epicyclic approximation. The 
current code, Epic5, derives the response of interstellar matter, orig- 
inally in circular orbits, to the gravitational potential in a matter of 
few seconds, which makes it an efficient code for surveying a large 
parameter space. In section [2] we present the analytical solution 
for the epicyclic theory with a generic potential, and we describe 
our code Epic5. We compare results generated by Epic5 with |Wada| 
(19941 in section [3] and in section [4] we demonstrate a standard 
case of a barred galaxy. In section [5] we compare results obtained 
by Epic5 with observations and simulations made for the galaxy 
NGC 1365. Finally, we conclude in section[6] 



2 METHOD 

The epicyclic description of nearly circular stellar orbits in a cir- 
cularly symmetric galaxy was developed by B. |Lindblad| {1927} 
in order to theoretically explain the observed velocity ellipsoid in 
the Milky Way galaxy. Later on, B. Lindblad (19581 introduced 
a rotating perturbing potential in the theory and pointed out the 
resonances that carry his name. P.O. Lindblad & P.A.B. Lindblad 
(1994), as well as Wada ( 19941, introduced gas dynamical friction 
in the epicyclic approximation. This form of damping the motions 
would make the first order approximation valid also over the Lind- 
blad resonances, but still not over the corotation resonance. 
The potential we are considering in our problem is: 



0(r,ff) = Oo(/) + 1>iO,#) 



(1) 



where (r, S) are the polar coordinates in a frame co-rotating with the 
potential at the pattern speed frequency Cl p . <t>o(r) is the axisymmet- 
ric potential, and 9) is an arbitrary perturbing potential devel- 
oped as a Fourier series 



<S>i(r,0) = -J^ [P c ,„,(r) cos (m 9) + P s , m (r) sin (m 9)] (2) 



or 



<S>(r, 9) = a> (r) - J] ¥„,(r) cos m (9 - & m (r)) (3) 



where 



= P;,„, + I*, and tan( ml ?,„) = P c ,,„/P s ,„, (4) 

We consider small deviations, £ and r], from circular motion 
following the usual notation, where 



r = r + £ 



1 



9 + (Q - Sl„) t + —r] 
ro 



(5) 
(6) 



and Cl(r ) is the angular frequency of circular motion. In addition, 
we introduce a frictional force proportional to the velocity devia- 
tion from circular motion with the coefficient -2A. As Wada ( 1994 1 
points out, this is analogous to the Stokes' formula (eg. |Yih|1979") 



p. 365) where the drag on a sphere moving slowly in a viscous fluid 
is proportional to the first power of the velocity. The equations of 
motion, as derived in appendix|A] can then be written as: 



I + 2A£, - 2Q.TJ - 4QAf : 



d<S>i 
'~d7 



= ^] [Cm cos m{9 - )?,„) + E m sin m{6 - §„)] 

m=l 

(7) 

i) + 2Q£ + 2Xi] + 4AA£ = - - ^ = - V D„, sin m(6 - & m ) (8) 
r ad t— 1 

where 

r = m - n = m — - E =rri¥ — 
ar r ar 

The solution for the motion of the guiding center will then be: 

n 

? = ^ [dm cos m(6 -& m ) + e m sin m(6 - !?„)] (9) 

m=l 
/t 

r] = ^ [gm smm(9 - & m ) + f,„ cosm(0 - # m )] (10) 
where the amplitudes d m , e m , g m and f m are given in appendix [A] 

2.1 The code: Epic5 



The code Epic5 is an extension of the code developed by Lindblad 
& Lindblad (1994) and computes the analytic solution and (10) , 
as well as its corresponding density and velocity maps, generated 
by the arbitrary potential. Epic5 derives the axisymmetric poten- 
tial, Oo, from a given rotational velocity curve and the perturbing 
potential, Oi, is introduced by its Fourier decomposition, P cm and 
Ps,„ (eq.0. 

In addition to the velocity curve and the parameters of the po- 
tential, another free input parameter is a constant pattern speed of 
the potential, Q p . The damping coefficient of the frictional force, 
A, which was introduced in the analysis to simulate orbits of in- 
terstellar matter, is assumed to have a linear tendency with radius. 
This linearity is determined by the given values of A at the outer 
inner Lindblad resonance (oILR) and at the outer Lindblad reso- 
nance (OLR) of the system (it is taken at kpc if no ILR is present 
and if no OLR is present, it is taken at r max in the rotation curve). 
These two values are also given as input parameters in Epic5. This 
definition of the damping coefficient, A, that allows it to vary along 
the radius, is desirable since it should depend on the gas density. 

As can be seen from eqs. |A13|A16| the introduction of the 
friction coefficient A has damped the amplitudes and eliminated the 
singularities at the Lindblad resonances where k 2 - (J^ = 0. How- 
ever, the singularity at corotation, (o m = 0, remains. As shown by 
|Binney &~T remaine (2008 Ch 3.3.3(b)), the stellar motions close 
to corotation turn into pendulum like oscillations around the La- 
grange equilibrium points L 4 and L 5 . With increasing distance from 
corotation, the angular amplitude increases until the orbit reaches 
the equilibrium points L\ and Li and flips into a circumcentral or- 
bit (see Figs. 16 and 26 in P.A.B. |Lindblad et al.||1996| hereafter 
LLA96). In the simulations of LLA96 the points L 4 and L 5 rep- 
resent density minima and L\ and Li density maxima along the 
corotation radius. 

Actually, our equations of motion j9} and \\Q\ are no longer 
valid around corotation because when inserting 9 from eq. |6]l into 
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Figure 1. Damping of the amplitudes f and r\ at corotation. Solid line: Vari- 
ation of the amplitudes with at according to the first order epicyclic theory 
with e = 0. Dashed line: Variation of the amplitudes introduced in Epic5. 



the differential equations l|7j and ([8} we have assumed rj/r « (Q.- 
£2 p ) which is not true close to corotation where SI - £l p approaches 
zero. We then see from eq. |8} that for small £ , £ and r\ we get a 
pendulum like oscillation in t], and we should not get a linearized 
equation. 

P.O. |Lindblad| ([T960 1 has pointed out that a time dependent 
potential can damp all resonances, including the corotation reso- 
nance. For a potential that varies as e yl , besides the twisting of the 
orbits around corotation, tu m in the denominator of eqs. ( |A13} to 



( A16 1 is replaced by V( w m + 7 2 )- m similarity to this, in Epic5 the 



corotation singularity is avoided by replacing the factor to,,, in the 
denominator of eqs. ( |A13^ - (AT6| by co € , where 



(11) 



where e„, = me and e is an additional input parameter in the code. 
As illustrated in Fig.[T| the amplitudes of the motion of the guiding 
center, growing large and of opposite sign on both sides of corota- 
tion, are smoothed and brought to zero at the exact resonance dis- 
tance. l/tJ £ reaches a maximum of l/2e„, at the distance a>,„ = e,„ 
from corotation. This means that the stellar elliptical orbits, trapped 
around corotation, in our case collapses to circular rotation at the 
corotation radius. This may not be an entirely inappropriate approx- 
imation to the solution in the corotation case. 

The derivation of the density is made applying the continuity 
equation. Using the equation (F-5) in Binney & Tremaine (20081, 
in our notation the ratio between the perturbed and unperturbed 
surface densities will be 



P_ 



~rd9 



8r 



(12) 



where we have assumed the unperturbed surface density p to have 
a flat distribution. 

Epic5 generates a number of supporting plots and diagrams 
like the rotation curve, the circular frequencies showing the reso- 
nances, to vs k, the perturbing potential, the total potential, radial 
and tangential perturbing forces, the phase variation of the perturb- 
ing potential, orbits, densities, as well as line of sight velocity fields 
for various position angles of the line of nodes. 

The resolution in r required by Epic5 in the tabulated values 



of the Fourier components of the perturbing potential is related to 
the increase of the phase shifts &„, with r. It is convenient to ex- 
press this relation in terms of the wavelength of the mode m = 2 
in the spiral structure. If A m is the wavelength of mode m, we have 
approximately for a tightly wound spiral 



dr 



In 
mX m 



n 



(13) 



Our requirement for Epic5 to interpret the potential tables P c , m (r), 
P s , m (r) is A(m& m ) < n/2 or 



Ar < 



2n 



(14) 



where Ar is the radial increment in the potential table and n is the 
maximum value of m. 

Finally, we can choose between counter-clock and clock wise 
rotation for an easier comparison with observed cases. 



3 Epic5 VS WADA'S (1994) MODEL 



|Wada| ( |1994| hereafter W94) independently provided an analyti- 
cal model, similar to Epic, representing the behavior of a non-self- 
gravitating gas in a rotating potential with a weak bar-like distor- 
tion. Like the original version of Epic ( Lindblad & Lindbla d]1994| > 
there was no solution for the corotation region. 

In the model Wada used the Toomre potential for the axisym- 
metric potential, and for the perturbation the potential given by 
|Sanders|(l977l > 



$>(R, if/) = *oW(l + cos 2t//) 



& <D b = e a 



R 2 



(R 2 +a 2 Y 



(15) 



(16) 



In these equations, v max is the maximum rotational velocity, 
a is the core radius and £o is a parameter which represents the 
strength of the bar potential. 

We test our code by using the same potential, circular fre- 
quency and parameters as W94. W94, however, chooses the friction 
to be proportional to the velocity in the radial direction only. 

|W94| considers a core radius, a, equal to 1 and a maximum 
rotational velocity, v raax , equal to (4/27) 1 ' 4 , simplifying the expres- 
sions of ® an d W94 assumed a weak bar strength, e , equal to 
0.05 and a pattern speed equal to 0.1. 

We use in this section the same potentials and derive the input 
needed for Epic5 from them (the rotational velocity curve and the 
Fourier components of the perturbed potential). We have also used 
the same pattern speed of 0.1, a friction coefficient of 0.02 at the 
oILR and 0.01 at the OLR and a corotation softening coefficient of 
0.02. We have located the bar horizontally, equally to its position 
in W94, and consider a counter-clockwise rotation. 

Using these input parameters for our code, we have derived 
the resonance radia at the same radia as W94: inner/outer Lindblad 
resonance (ilLR and oILR) at 0.90 and 2.22 respectively, corotation 
(CR) at 4.53 and outer Lindblad resonance (OLR) at 6.09 (see Fig. 
1 in W94). We present the orbits generated by Epic5 in the upper 
panel of Fig. [2] corresponding to the ones presented in Fig. 4 in 
W94. The leading arms around the ilLR and trailing arms around 
oILR are well presented in both cases. Epic5 damps as well the 
orbits around corotation, allowing us to present orbits until radia 
further than CR (see middle panel in Fig. [2j. 
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Figure 2. Results from Epic5 for the W94 case. Upper graph: Orbits gener- 
ated by Epic5 inside corotation. Middle graph: Orbits generated by Epic5 in- 
cluding the outer Lindblad resonance (OLR). Lower graph: Density con- 
tours given by Epic5 according to eq. \\-\ Contours at 1.01, 1.05, 1.1, 1.2, 
1.3 and 1.4. Red circles show the positions of resonances ilLR, oILR, CR 
and OLR. The rotation is counter clockwise. 
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Figure 3. Rotation curve given by two superimposed exponential disks with 
masses of 11.5 X 10' M and 2.87 X 10 9 M and with scale lengths of 
0.194 kpc and 1.15 kpc respectively. 

Epic5 estimates the density map generated by the input poten- 
tial in the lower panel in Fig. [2] showing the continuation of the 
spiral arms out to the OLR. 



4 Epic5: A DEMONSTRATION 

We demonstrate in this section a case to illustrate what Epic5 is 
able to produce. The rotation curve in this case, shown in Fig. [3] is 
that given by two superimposed exponential disks with masses of 
1 1.5 x 10 9 M G and 2.87 x 10 9 M Q and with scale lengths of 0. 194 kpc 
and 1.15 kpc respectively (reasonable values for the inner part of 
barred galaxies, Lin dblad et al.|2010| >. To this, is added an ad hoc 
perturbing bar potential, described by its Fourier components of 
cosine of 2, 4, 6, 8 and 10 times 6 for each value of the radius, as 
seen in Fig. [4] (left graph). This results in the total potential shown 
in Fig. [4] (right graph). The relative bar strength IdhUn/IOolmin is 
equal to 0.026. 

Fig. [5] shows the variations of £2, CI ± k/2 and CI ± k/4 with 
radius. Our assumed pattern velocity of 70 km s kpc -1 , which 
places the oILR close to the minimum of the perturbing potential, 
is shown as a straight line in the figure, and the resonances are 
shown by the vertical lines. The ilLR lies very close to the center, 
oILR at 1.33 kpc, CR at 2.26 kpc, and OLR at 3.03 kpc. 

Computing the orbits, we have assumed the damping coeffi- 
cients constant at A = 60 km s kpc~' and e = 60 km s -1 kpc -1 . 
The resulting orbits are seen in Fig. [6] The damping coefficients 
have been adjusted such that the orbits, presenting a laminar flow 
pattern, do not cross. This is also a condition for Epic5 to be able to 
produce meaningful density and velocity maps. This choice means 
that these coefficients and the strength of the bar are correlated. 
The densities, as given by Epic5 are shown in Fig. [7] (left graph), 
and the line of sight velocities in the plane of the galaxy as seen 
with a position angle of the line of nodes of 0° is given in Fig. [7] 
(right graph). 

We see how the orbits well inside the oILR are orientated per- 
pendicular to the bar but twisting against the direction of rotation. 
At the oILR they are tilted about 45° against the position of the 
bar, and close to corotation they are elongated along the bar. The 
orbits show sharp kinks on the leading side of the bar. This is con- 
nected to large velocity jumps at the corresponding positions and 
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Figure 5. Angular velocity, CI, against radius derived from the rotation 
curve (solid line). The dashed line shows n - k/2, the dotted-dashed line 
shows fl + k/2, the dotted and dashed-triple dotted lines show CI ± k/4. The 
pattern speed of 70 km s _1 kpc~' is shown by the horizontal solid line, and 
the resonances are indicated by the vertical lines. 

strong enhancements of the intensities parallel to the bar on the 
leading side from the ILR to CR. Such dust lanes along the leading 
side of the bar are well-known common features in barred galax- 
ies (e.g. Athanassoula 19921. From corotation the response splits 
and a trailing spiral arm continues out to OLR. The very faintest 
density contours close an ellipse between CR and OLR elongated 
perpendicular to the bar. 

Using the same case, we have varied O,, to observe changes in 
the response. The density and velocity maps generated by Epic5 for 
Q. p equal to 60 km s _1 kpc -1 and to 80 km s _I kpc -1 are shown in 
Fig. [8] We have used a constant damping coefficient A and e equal 
to 60 km s kpc~' for the slow and fast potential. The trailing arms 
extending from CR and going towards OLR are shorter and fainter 
as the potential gets slower. Moreover, a faster bar generates more 



ORBITS 

-2 2 




-2 2 

Figure 6. Orbits generated with Epic5. The red circles mark the oILR, coro- 
tation and OLR. The rotation is clockwise. 

prominent, and also more tightly wound, dust lanes inside CR, with 
larger velocity jumps, as shown in Fig. [8] right graphs. When the 
pattern speed gets larger, the resonance distances decrease. As seen 
from Fig. [4] (left graph), the CR and OLR move to stronger per- 
turbing amplitudes, which is partly the reason for the change of 
response. With an observed bar potential and observed structures 
similar experiments should help to estimate the pattern velocity and 
the true positions of the resonances. 
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Figure 7. Left graph: Density contour p/pi in the plane of the galaxy. Contours at 0.99, 1.01, 1.02, 1.05, 1.1, 1.2, 1.4 and 1.7. Right graph: Velocity field in 
the plane of the galaxy with a position angle of 0° and an increment of velocity between contours of 50 km s~ 1 . The direction of rotation is clockwise. 



5 Epic5 APPLIED TO NGC 1365 

To evaluate the limitations and abilities of Epic5 we want to com- 
pare it to hydrodynamic simulations of a nearby barred galaxy. For 
this, the strongly barred Seyfert galaxy NGC 1365 (Fig. [9} would 
be suitable. This galaxy has been extensively observed (see the re- 
view by P.O. Lindblad 1999) and a hydrodynamic simulation of 
the gas flow pattern was derived by LLA96. Also Zan mar Sanchez| 
et al. (2008 ) observed and modeled the gas flow in the bar. 



A total hybrid velocity field has been derived by Lindblad et al 



1996[l based on the HI velocity field of Jorsater & van Moorsel 



1995[l and complemented with optical long slit spectra. We have 



here accepted the rotation curve based on this as derived by LLA96 
(Fig.|10[>. It ag rees closely with the rotation curve given by |Zanmar| 
Sanchez et al. ( 2008 ). We have derived the perturbed potential from 
the total surface density following the analysis using Bessel func- 
tions described in |Binney & Tremaine|(2008| Ch. 2.6.2) and LLA96 
(their Appendix A). For this analysis, we used the Fourier series de- 
composition of the total surface density, obtained by LLA96 from 
a J-band image as given in their Figs. 7 and 9, and a normalized 
triangle density distribution along the vertical direction, with a zo 
scale parameter equal to 1 kpc. We smoothed the rotational velocity 
curve and the Fourier components along the radius and interpolated 
the curves in order to achieve enough resolution for the A 2 criterion 
of the spiral potential (eq.|14|). We use a clock-wise rotation in order 
to reproduce the observed galaxy. 

When comparing with observations, we have taken the results 
given by Epic5, located in the galactic plane, and transformed them 
into the sky plane. For this purpose, we have used a position an- 
gle, PA, equal to 40" and an inclination angle, i, of 42" (Jorsater 
& van Moorsel 1995, also used in LLA96). In addition, we have 
assumed the systemic velocity to be 1632 km s~' and a scale of 
0.1 kpc/arcsec. 



5.1 Comparing with observations and simulations 

The hydrodynamic simulations presented in LLA96 are performed 
using a flux-splitting 2nd-order code FS2, solving the flow of an 
ideal isothermal non-viscous gas ( [van Alba da 1985). LLA96 ap- 
plies two different potentials derived from NGC 1365: bar only and 
bar + spiral. We have followed the same procedure in order to com- 
pare the results. 



5.1.1 Bar potential 

When deriving the bar only potential, only the cosine contributions 
to the Fourier series has been considered, and only up to a radius 
of around 100 arcsec, where the phase shift of the different modes 
start a steep increase (see Figs. 8, 9 and 1 1 in LLA96). The bar is 
positioned horizontally along the x-axis (Fig.[TT| upper left). 

As in LLA96, we have used a pattern speed of 20 
km s~ 



kpc . As illustrated in Fig. 11 (lower left), this pattern 
speed fixed the resonance radia to: the inner inner Lindblad reso- 
nance (ilLR) at 0.4 kpc, the outer inner Lindblad resonance (oILR) 
at 2.7 kpc, the corotation radius (CR) at 14.4 kpc and the outer 
Lindblad resonance (OLR) at 21 kpc. We have used a constant fric- 
tion coefficient A of 7 km s -1 kpc -1 , and a corotation damping co- 
efficient £ of 5 km s -1 kpc -1 . We are using an amplitude for the 
perturbing potential A bal = 0.5, half the amplitude used in LLA96. 
Thus, we have a potential with a bar strength, |®i| m i n /|®olmin, °f 
0.01. Increasing the bar amplitude, means that we have to increase 
the damping coefficients correspondingly to avoid self-crossing or- 
bits, that are not realistic in a gaseous flow, and the results will be 
nearly the same. This may imply that the bar in NGC 1365 is actu- 
ally too strong to be dealt with Epic5. 

Using these input parameters, we obtained with Epic5 the or- 
bits, densities and velocity field shown in Fig.[TT]and Fig. [12] The 
orbits are seen face-on. The density and velocity maps are projected 
on the sky plane, for comparison with the results of the simulations 
derived in LLA96 as well as with the observations. 
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Figure 8. Upper and lower left graph: Density contour p/pi in the plane of the galaxy, for Ci p equal to 60 km s kpc~'and 80 km s _1 kpc -1 , respectively. 
Contours at 0.99, 1.01, 1.02, 1.05, 1.1, 1.2, 1.4 and 1.7. Upper and lower right graph: Velocity field in the plane of the galaxy with aposition angle of 0° and 
an increment of velocity between contours of 50 km s . Upper graph for a n ; , equal to 60 km s kpc~' and lower for 80 km s kpc - '. The direction of 
rotation is clockwise. 



We see from Fig. [TT| (lower right) that the orbits in the re- 
gion around oILR take elliptical shape, twisting in the counter ro- 
tation direction from being orientated at large angles to the bar, 
over 45 degrees at the oILR, to more elongated with the bar outside 
this resonance. Inside corotation matter circulates clockwise with 
respect to the bar, and the twisting and increasing deviation from 
elliptical shape gives rise to density increase forming a small trail- 
ing spiral across the oILR. This straightens out to a lane along the 
leading edge of the bar with sharp shock-like velocity jumps over 
the lane. 

The influence of the inner 1:4 resonance, that occurs around 



R = 9 kpc, can also be seen in the orbits. Even the 1:6 resonance 
further out is discernible in the slight crowding of orbits. Strong 
crowding of orbits are seen at CR at the Lagrangian points LI and 
L2 at the end of the bar and, with slightly spiral form, at the OLR 
in a direction perpendicular to the bar. 

Compared to the orbits of the BM model of LLA96 (their 
Fig. 16), we see that between oILR and CR the BM orbits are more 
elongated and more orientated along the bar than ours. The 1:4 res- 
onance is clearly seen also here. Around corotation the BM model 
displays the pendulum like motion referred to before and which is 
not reproduced by the linear approximation of Epic5. 
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Figure 9. Optical image of NGC 1365 obtained at the prime focus of the ESO 3.6 m telescope. 



The interstellar matter density map (Fig. [12] left graph), now 
projected on the sky plane, should be compared with the optical 
image of the galaxy (Fig.|9| and LLA96 (their Fig. 15). As expected 
from the orbits, the map shows a lane on the leading side of the bar 
that bends into spiral shape over the oILR. At the position of the 1 :4 
resonance a faint spur is seen, a strong arm appears at the end of the 
bar along the corotation radius, and a faint trailing arm extends out 
to the OLR. All these features have their correspondence in the BM 
density map of LLA96 (their Fig. 15). Fig.[T2] right graph, shows 
the innermost region of the EPIC5 density map, where, as expected, 
a small nuclear leading spiral is seen across the ilLR. 

The velocity map, to be compared with LLA96 (their Fig. 20) 
is shown in Fig. [12] (middle graph). The map is projected on the 
sky plane with a line of nodes in PA 40° and an inclination between 
the galaxy symmetry plane and the sky plane of i = 42°. The main 



difference, comparing the Epic5 result against the observations and 
simulations, is that the twisting of the zero velocity contour towards 
the bar is much weaker in the Epic5 case, which is connected with 
the different shapes of the orbits in the bar discussed above. This 
may partly be due to the difficulty of Epic5 to accept a very strong 
bar, as in NGC 1365, without getting crossing orbits. The behav- 
ior of the velocity contours inside the oILR are similar to that in 
LLA96, and the sharp velocity jumps, of the order of 100 km s , 
in the lanes on the leading side of the bar are seen in all cases. The 
twisting of the contours along the spiral arms with velocity jumps 
on the inner side of the arms inside corotation are seen both in the 
simulations and in the Epic5 results. 
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Figure 10. Rotational velocity curve of NGC 1365 (from LLA96). 
5.1.2 Bar + Spiral potential 

In the case of the bar + spiral potential we have taken the sine 
and cosine contribution of the Fourier series of the first three even 
modes until a radius of around 220 arcsec, just as described in 
LLA96. The resulting potential is seen in Fig. [T3] We have used 
a pattern speed of 18 km s -1 kpc -1 (found for this potential by 
LLA96) and placed the bar horizontally along the x-axis. This pat- 
tern speed makes the resonance radia to be: the inner inner Lindblad 
resonance (ilLR) at 0.38 kpc, the outer inner Lindblad resonance 
(oILR) at 3.1 kpc, the corotation radius (CR) at 15.8 kpc and the 
outer Lindblad resonance (OLR) at 23.5 kpc. We have chosen the 
coefficients A and e and parameter A], m the same as in the bar only 
case. The main difference between this case and the bar only case 
is that the spiral features are led much easier through corotation, 
which can be seen both in the orbits, Fig.^Jlower right graph, and 
the densities, Fig. [14] upper left graph. This was noticed also by 
LLA96. 

Fabry-Perot interferometry of NGC 1365, covering the Ha 
regime has been obtained by Zanm ar" Sanchez et al.| {2008}. By 
kind permission we have got access to their total Ha intensity data, 
which are reproduced in Fig. |15| (left map). In Fig. |15| (right map) 
we have for comparison color coded the Epic5 intensities from Fig. 
[T4| Clearly the Epic5 model agrees well with the Ha structure. The 
hot spots in the nuclear region and the active Seyfert nucleus is of 
course not reproduced by Epic5. 

The velocity field obtained with Epic5 is shown in Fig. [14] 
together with the optical radial velocity field of NGC 1365 (in- 
ner 100 arcsec) taken from LLA96. In the innermost part, well 
inside the oILR, the contours indicate rapid rotation in nearly cir- 
cular orbits. Between oILR and CR we again see how the isove- 
locity contours are concentrated to the leading edge of the bar. In 
the Epic5 case the contours are squeezed to a shock-like jump of 
about 100 km/s. The optical map cannot make such a sharp jump 
because it is smoothed due to the limited spatial resolution given 
by the limited set of long slits on which it is based. Over corotation 
we see similar shock behavior on the inner side of the spiral arms. 
The residual map obtained by subtracting the rotational velocity, 
Fig. [TO] from the Epic5 model velocity field is presented in Fig.|16| 



In this Fig. 16 is also presented the residual velocities taken from 
|Jorsater &~v an Moorsel ( 1995 1. They subtracted from the observed 
HI velocity map, the rotational velocity derived from their obser- 
vations considering a warp disc model after a radius of 250 arcsec. 



The general behavior is similar in the two cases with shocks along 
the bar and spiral arms. It is informative to compare the Epic5 map 
with the orbital map in Fig.[T3](lower right) to envision how the or- 
bital circulation results in the observed line of sight velocity map. 



6 CONCLUSIONS 

Epic5 cannot fully replace elaborate hydrodynamical simulations. 
Epic5 solves the equations of motion in a solid body rotating, 
time independent galactic potential. This is done in the first or- 
der epicyclic approximation. Epic5 computes the deformation of 
the initially circular orbits of the guiding center, the velocity fields, 
and the structure created. Epic5 pictures a steady state created by 
a weak perturbing potential, it does not give the time evolution of 
galactic structure. Due to the degeneracy of the strength of the grav- 
itational potential and damping coefficients, we can only derive rel- 
ative values of the density. 

Our comparison with a more elaborate simulation, as well as 
observations, show that it still can reproduce structural and dy- 
namic properties. One of the advantages with Epic5, involving an 
analytic solution, is the speed by which it can survey the parameter 
space. Further, the possibilities, by varying parameters, to isolate 
artificially various agents important for forming structure of galax- 
ies will be of importance for understanding the formation of galac- 
tic structure. Epic5 will be a useful complement to time consuming, 
elaborate simulations. 

We would like to thank to J. A. Sellwood and R. Zanmar 
Sanchez for letting us use their Ha data of NGC 1365. We also 
want to thank to the reviewer, Gene Byrd, for his useful comments. 
NP-F acknowledges financial support from NOTSA, and KF is sup- 
ported by the Swedish Research Council (Vetenskapsradet). 
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Figure 12. Left graph: Density ratio, p/po- Contours at 1.06, 1.2, 1.4 and 1.6. Middle graph: Line of sight velocity map derived by Epic5 using the position 
angle and the inclination of NGC 1365. Contours every 30 km s~', and thick solid line for the systematic velocity of the galaxy, 1632 km s -1 . Right graph: 
Density ratio, p/po, for the inner 4 kpc. Contours at 1.02, 1.04, 1.06, 1.2, 1.4 and 1.6. Maps seen projected on the sky plane. 
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Figure 13. Upper left graph: Bar + spiral perturbing potential contours of NGC 1365. Upper rigth graph: Total potential contours of NGC 1365 corresponding 
to a perturbing strength of the potential, |Oi| m i n /|Oo| m i n , equal to 0.01. Red curves: inner and outer Inner Lindblad Resonance (ILR), corotation and Outer 
Lindblad Resonance (OLR). Lower left graph: Phase shift of the perturbed potential, &„,(r), presented in the appendix^] for the three considered modes of its 
Fourier development. Thick solid line correspond to m equal 2. Lower right graph: Orbits seen as face-on generated by Epic5 using the bar + spiral potential 
and a constant friction coefficient, A, equal to 7 km s - ' kpc~' and a corotation damping coefficient, e, of 5 km s~' kpc~' . 
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Figure 14. Upper left graph: Density ratio, p/po, for bar + spiral potential. Contours at 1.06, 1.2, 1.4 and 1.6. Upper right graph: Line of sight velocity map 
derived by Epic5 using the position angle and the inclination of NGC 1365. Contours every 30 km s , and thick solid line for the systematic velocity of the 
galaxy, 1632 km s . Lower left graph: Optical radial velocity field in the inner 100 arcsec of NGC 1365 taken from LLA96. Lower right graph: Line of sight 
velocity map derived by Epic5 for a bar + spiral potential in the inner 100 arcsec. Maps seen projected on the sky plane. 
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Figure 15. Upper left map: Ha intensity map observed using a Fabry-Perot interferometer by Zanmar Sanchez et al. (2008). Shown in arbitrary units. Upper 
right map: Model density ratio p/po obtained by Epic5 , eq. [12) . Ratios lower than 0.9 are in white color. 




Figure 16. Left graph: Velocity model generated by Epic5 for a bar + spiral potential minus the rotational velocity from Fig.|10| Den sity contour over-plot in 
solid black lines. Right graph: HI velocity field observations minus rotational velocity model, figure from Jorsater & van Moorsel ( 1995 1. The contours are 
from an optical image of NGC 1365. The NW side of the galaxy is the near one, and the rotation is clockwise, x-axis and v-axis in both graphs in arcsec. 
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APPENDIX A: PERTURBATIONS IN THE EPICYCLIC APPROXIMATION 

We present in this paper a solution of the motion of a mass particle in an axisymmetric potential, O (r), with a weak non-axisymmetric 
perturbation, Oi(r, 9). We assume that the potential is rotating with a pattern speed Cl p , and we study the system in a corotating frame. We 
express the potential in polar coordinates: 



<P(r, 0) = 4>o0) + 9) = 4> (r) - V <P m (r) cos m(0 - K(r)) (Al) 



m= I 



where &„, is the spiral phase. The rotation of the axisymmetric potential is about the center of gravity, which we have assumed to be at rest. 
Thus one has to be cautious in the presence of strong m = 1 perturbation, when the center of mass could be significantly displaced. 
In a corotating coordinate system we can write the equations of motions (see also Binney & Tremaine 2008 p. 189): 



<9® 
'Ik 
1 dO 



+ 2rOCl p + rCl l p (A2) 



r9 + 2r9 = 2rQ„ (A3) 

r 89 p 



We introduce f (f) and 77(f) as deviations from circular motion and write: 



r = r +( (A4) 
1 

= 6/ + (CI - O p ) t + — 77 (A5) 
r 

where CI is the circular angular velocity at a radius r . Assuming £ and 77 to be small, and linearizing the equations of motion by neglecting 
higher order terms of £ and 77, we get: 

'£ - 2Q.T] - 4QA£ = -^P- (A6) 
or 

i7 + 2£^ = --^ (A7) 
r aG 

where A is the Oort constant, A = -r/2 dSl/dr. 

To describe the motion of a gaseous medium, we introduce a frictional force proportional to the deviation from circular motion with a 
damping coefficient, 2A (see Wada 1994, Lindblad & Lindblad 19941. With the introduction of this frictional force the equations of motion, 
eq. \A6\ and |A7j, are transformed into the equations ( |A8[ l. 

f + 2^ - 2fl77 - 4QA£ = -—! = > [C m cos m(0 - &„,) + E m sin m(6 - &,„)] 
or -<H 

1 ao, A 

77 + 2Q£ + 2Af] + AXA£, = = - > D,„ sin m(0 - & m ) (A8) 

r 89 

117=1 

where 

r = ZlE- n = m — - E =nW — 
ar r ar 

Provided that we are not close to corotation, where CI = Q. p , we can assume that 77/ro << (CI - Sip), what makes 

e~e + (Ci-ci„)t (A9) 



with k = 4C1 - 4£X4, we can write the full solution of the now linearized eqs. 1 A8 1 as 



£ = ce~ aM C0SK(t - tp) + ^ [d,„ cosm(6 - #,„) + e„, s'mm(9 - & m )] (A10) 

m=l 

2CI " 

77 = ce^ A ' su\K(t - to) + / [g,„ sinm(6 - 1?,,,) + f,„ cosm(9 - #,„)] (All) 

a and /3 are functions of CI and k, and c is an arbitrary constant. The second terms on the left side of eqs. {A10| and |ATTJ give the forced 
oscillation due to the perturbing force. The first terms give the damped oscillation with the epicyclic frequency k around these guiding centra. 
We will leave out these latter terms in what follows and just consider the motions of the guiding center. 
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We introduce these eq. jA10| l and eq. {ATT} into |A8j, getting a system of equations with the 4n unknown amplitudes, d m , e m , g,„ and 

fm • 



-(to 2 m + 4CiA)d m - 2co m Clg m 
-(col, + 4nA)e,„ + 2co m Clf m 
2co m Cld m + co 2 m g m 
2co,„Cle m - coljm 
where co„, = m (CI - Cl p ). The final solution is shown below. 



co s [(k 2 -coI) 2 + M 2 (k 2 +co 2 ,) + 16A*] 



Cm 2AC0 ln C m 

= E m + 2Aco„,d„, 

= D m - 2Aco,„f„, + 4AAe m 

= -2Aco m g„, - 4AAd m (A 12) 



(A 13) 



(k 2 - (4)(to m C m + 2C1D,„) - 2A(k 2 + to 2 m )E m - 4A 2 (co m C m - 2ClD m ) - SA 3 E, 



-(k 2 - col) [2« {to m C m + 2CIDJ - (k 2 - co 2 m ) D m \ + 4A [{k 2 - co 2 m ) A + 2co 2 m £l] E m 

" col[(K 2 -co 2 m ) 2 + M 2 (K 2 + col) + l6A*} 

-4A 2 \(co m C m + 2ClD m ) (2C1 - 4A) - (co 2 , + 3k 2 - 8fl 2 ) D,„] + \6A 3 AE m 
_i t i ' 1 (A14) 

co 2 e [{k 2 - co 2 J + M 2 {k 2 + co 2 ,) + 16A*] 



col {« 2 - <4) E m + 2A [{k 2 + col) (oo m C m + 2ClD m ) - 2Ci (k 2 - col) D m ] 


j ~ 4/l 2 ai 2 ,£ m + 8A 3 co m C m 


col[(K 2 -colf + M 2 {K 2 + col) + l6A*) 


1 



(A 15) 



2co m [k 2 - co 2 ^j C1E„, + 2A ^4co m Cl (co m C,„ + 2ClD m ) + (k 2 


- co 2 m ) (2AC m - co m D m )] + U 2 co m E m (CI - 2A) + U 3 (2AC m + co m D m ) 




) 2 + U 2 (k 2 + col) + 16-* 4 ] 



/» = ' F 5 — : ' S ( A16 ) 

coI[(k 2 -coI) 2 + 8A 2 (k 2 + coI) + 16A<] 

At the Lindblad resonances k 2 - col = 0> anc ' we can see now me friction terms are damping the resonances when A + 0. At corotation 
co„, = 0. To damp this resonance, as explained in Section[2]in the paper, we have replaced co„, in the denominator of eqs. ( |A13} to l |A16[ l by 
co E , where 

1 co.„ 

— = "rT^ (A17) 

where e,„ = me and e is an additional parameter. 
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